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Abstract 



We calculate the production oi a W boson and a single b jet to next-to-leading 
order in QCD at the Fermilab Tevatron and the CERN Large Hadron Collider. Both 
exclusive and inclusive cross sections are presented. We separately consider the cross 
section for jets containing a single b quark and jets containing a bb pair. There are 
a wide variety of processes that contribute, and it is necessary to include them all in 
order to have a complete description at both colliders. 
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1 Introduction 



Many signals of new physics at hadron colliders involve a weak vector boson {W, Z) plus jets 
containing heavy quarks (c, fe). For example, the top quark was discovered via the signal 
W + ij, with at least one 6 jet [H [2]- More recently, evidence for single top production has 
been presented via the signal W + 2j, with at least one 6 jet [3], HI [5l [6] . The Higgs boson 
could manifest itself via the same signal, from the production process qq — > Wh, followed by 
h bb [3, El E]. The Higgs boson could also appear in the signal Z + 2j with at least one b 
jet, via qq —>■ Zh [10, ITT] - 

Most calculations of the background processes that give rise to W, Z + nj {n = 1,2), 
where one or more jets contain heavy quarks [Q = c,b), have been completed at next-to- 
leading order (NLO) [121 [13]. There exist NLO calculations of [H], ZQQ [151 [I6], ZQj 
[n], M/'c [IHIIE], Wbb [201 EH, and ^bj [22]. An obvious omission from this hst is Wb, that 
is, W + Ij with at least one b jet at NLO. It is the goal of this paper to fill in that gap. 
This NLO calculation can be used to normalize the cross section from a leading-order event 
generator such as ALPGEN [23] or MadEvent [21]. 

It may seem surprising that the NLO calculation of Wb has not already been done, given 
all the other NLO calculations listed above. The reason for this, as we will discuss shortly, 
is that it is essential to do this calculation with a finite 6-quark mass. In contrast, most of 
the above-mentioned calculations were done with a vanishing heavy-quark mass, with the 
justification that the quark mass is negligible at high transverse momentum (pr)- The ability 
to do NLO calculations with a finite heavy-quark mass for this class of processes was only 
recently demonstrated, for Wbb, in Ref. [21] (and for Zbb in Ref. [IE]). We will use this 
calculation, together with the NLO calculation of Wbj [22], to generate the NLO calculation 
of Wb, including the effect of the 6-quark mass. 

We discuss the details of the calculation in Section [2l We then present results in Section 
|3] and conclusions in Section HI 

2 Wb at NLO 

The leading-order (LO) processes for the production oi a. W boson and one jet containing 
a b quark are shown in Fig. [1] In both cases there are two partons in the final state, 
but we require that only one of them (which contains a b quark) reside at high transverse 
momentum, with pxj > 15 GeV at the Fermilab Tevatron {^/s = 1.96 TeV pp) and pxj > 25 
GeV at the CERN Large Hadron Collider (LHC, ^/s = 14 TeV pp). We also require this 
parton to lie within a pseudorapidity of \rij\ < 2 at the Tevatron and |?7j| < 2.5 at the LHC. 
Furthermore, we demand that two partons be separated by |Ai?| > 0.7; if they are not, 
then their four-momenta are added and they are considered as occupying a single jet with 
this four-momentum. These requirements are made to crudely simulate the acceptance and 
resolution of the detectors. They are listed in Table [H along with the parameters used in the 
calculations. In all calculations, the light quarks are summed over q = u,d,s,c, including 
CKM mixing. 

The LO process shown in Fig. [l](b) contains a b quark in the initial state, and requires 
further discussion [261 Elj- To understand the usefulness of this approach, consider the 
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Table 1: Cuts used to simulate the acceptance and resolution of the detectors, and parameters 
used in the calculations. 



Tevatron: p^j > 15 GeV Irjjl < 2 
LHC: pTj > 25 GeV \rij\ < 2.5 
\AR^l\ > 0.7 \ARbj\ > 0.7 



Mw = 80.44 GeV 
LO: CTEQ6L1 

Hf = Mw 
a^°{Mz) = 0.130 
a^^°{Mz) = 0.118 

= 8M^Gf/V2 
Vud = = 0.975 



mb = 4.7 GeV 

NLO: CTEQ6M [25] 

fiR = Mw 

a^^{Mw) = 0.132 

a^^O(Miy) = 0.120 

Gf = 1.16639 X 10-^ GeV-2 

Vus = = 0.222 



q 




Figure 1: Leading-order processes for the production of a W boson and one jet, which 
contains a b quark. 

alternative approach, shown in Fig. [2] [2H]- In that approach the fe-quark mass is kept 
nonzero, so one obtains a finite result even when one of the b quarks is emitted at zero 
Pf- However, although the collinear singularity is regulated by the b mass, one obtains an 
enhancement factor of In Mw/mb- Another power of this factor appears at each order of 
perturbation theory, degrading the convergence of the series. To ameliorate this, one sums 
this enhancement factor to all orders into a b distribution function, and uses this function 
in the LO calculation of Fig. The other big advantage of this approach is that the LO 

process of Fig. Wljo) is simpler than that of Fig. [2|, and hence a NLO calculation becomes 
tractable. This effectively allows to include a set of higher order corrections to the process 
of Fig. [2] that would appear only at NNLO in the fixed order calculation and will probably 
not be available for quite some time. 

There are a variety of processes that must be calculated at NLO: 

1. qq' ^ Wbb at tree level [Fig. [11(a)] and one loop (mj, 7^ 0) 

2. qq' ^ Wbbg at tree level (m^ 7^ 0) 
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Figure 2: An alternative way of calculating the process in Fig. [T](b). 

3. hq Wbq' at tree level [Fig. [11(b)] and one loop {nif, = 0) 

4. bq ^ Wbq'g at tree level (mb = 0) 

5. bg ^ Wbq'q at tree level {rrih = 0) 

6. gq Wbbq' at tree level (Fig. [2]) (m^ ^ 0) 

Processes 1 and 2 are calculated with a non-zero 6-quark mass, using the code developed in 
Ref. [21]. Processes 3-5, which involve an initial-state b quark, are calculated with = 0, 
using the code developed in Ref. [22]. Process 6, calculated with a non-zero fe-quark mass, 
can be taken from either code. We used this process to cross-check the two codes. We notice 
that the counterterm that subtracts from Process 6 the logarithmic terms already included 
in Process 3 has been added, in our calculation, to Processes 3-5, since it shares the same 
final state phase space configuration. 

In the formalism that our calculation is based upon [261 121], the 6-quark mass is set to 
zero in any process in which the b quark appears in the initial state (Processes 3-5). This is 
not a limitation of the formalism, nor is it an approximation. In all other processes, where 
the b quark appears only in the final state, the 6-quark mass is kept nonzero, although it 
is often a good approximation to set it to zero. For the calculation we are performing, it is 
essential to keep the 6-quark mass nonzero. This is because we demand that only one b jet be 
at high pti and we do not restrict the other b quark. Thus there are regions of phase space 
where the bb invariant mass is not much greater than 2mfe, in which case it is very inaccurate 
to neglect the 6-quark mass [22]. This issue arises already at leading order, as evidenced 
by Fig. [31 where we show the bb invariant mass from Process 1 for the Wb exclusive cross 
section at the Tevatron for both mi, = and = 4.7 GeV (with ARjj > 0.7). We see that 
the cross section (the area under the curve) is very sensitive to the 6-quark mass, even with 
the required separation of the two b quarks. 

All NLO calculations are done in the MS scheme. When we refer to the NLO Processes 
2, 4-6, it is understood that initial-state coUinear singularities are subtracted in this scheme. 
For massless quarks, the coUinear singularities are regulated dimensionally, while for massive 
b quarks (Processes 2 and 6) they are regulated using a finite b mass [22l [261 [2Z]- The 
calculation of Processes 3-6 was done using the Monte Carlo code MCFM [29| [30] . 
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Figure 3: The differential cross section for W plus one 6 jet {ptj > 15 GeV, |?7j| < 2), vs. the 
invariant mass of the b quark and the other b quark (outside the fiducial region and separated 
by A.Rjj > 0.7), at the Fermilab Tevatron {\^ =1.96 TeV pp). Only the contribution from 
the LO subprocess qq Wbb is shown. The (solid, blue) curve includes the b quark mass, 
while the (dashed, red) curve does not. 



3 Results 

In Table Ej we give the exclusive cross section for + 1 jet, where the jet contains a b 
quark, and there are no other jets present within the acceptance (listed in Tabled]). We list 
the cross section for the case where there is only one b in the jet (denoted Wb), and when 
there are two 6's in the jet (denoted W{bb)). The tagging probability for a jet with two 6's 
differs from that of a jet with one b [3l]. We notice that in the W{bb) case the two b quarks 
can be coUinear and give origin to large logarithms of the form \n{p^/mb). Our results do 
not contain a resummation of these logarithms and are therefore subject to some degree of 
uncertainty. We could have rejected the W{bb) configuration in our W + 6-jet NLO results 
or used a different jet algorith as suggested, for instance, in Ref. [32]. We prefer to keep 
it and quote it separately in order to make it available to different kinds of experimental 
analyses. Table [2] gives both the LO cross section (in square brackets)0 and the NLO cross 
section. The first number given is from Processes 1-2; the second from Processes 3-6; and 

^The contribution from Process 1 to the Wb cross section at the Tevatron corresponds, at LO, to the 
area under the (sohd, blue) curve in Fig. [3l 
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the third is their sum. In the case of Wb we also quote, in parenthesis, the contribution of 
Process 6 by itself. Indeed Process 6 is part of both the fixed order (1 + 2 + 6) and resummed 
(3 + 4 + 5 + 6) calculation^ and results could be quoted grouping Processes either way. 
In order to assess the impact of resumming initial state coUinear logarithms, however, one 
needs to compare the sum of Processes 1,2, and 6 to the sum of process 3, 4, and 5. This can 
be easily deduced from our Tables knowing Process 6 independently. In the case of W{bb), 
only Process 1 contributes at LO, and only Processes 1-2,6 at NLO. Thus we do not need to 
quote Process 6 separately. The NLO results correspond to a pure fixed order calculation. 
We see that at the Tevatron, Processes 1-2 make a much larger contribution to Wb than 
Processes 3-6, while at the LHC the two sets of processes make comparable contributions. 
This mirrors the results for the Wbj final state [22], first noted in Ref. [28]. The fixed 
order calculation (Processes 1 + 2 + 6) dominates at both the Tevatron and the LHC, but 
the corrections included in Processes 3+4+5 are as much as 25% at the LHC and therefore 
very relevant. We also notice that for both Wb and W{bb) Processes 1-2 dominate at the 
Tevatron, while at the LHC Process 6 is nearly half the size of Processes 1-2. Therefore it is 
very important to include higher order corrections to this process, as it is achieved by adding 
to the fixed order result Process 3, 4 and 5 (see discussion in Sec. [2]). Finally, we notice that 
for W{bb), the NLO cross section is significantly larger than the LO cross section, while the 
NLO correction is modest for Wb. 

Also given in Table [2] is the exclusive cross section for Wj (j = u,d, s,c, g), both at LO 
and NLO [33l [M] . These numbers are useful to compute the fraction of + Ij events in 
which the jet contains a b quark. This fraction is around 0.7% at the Tevatron and 0.8% at 
the LHC. 

We give in Table [3] the inclusive cross sections for W + Ij + X, where the jet contains a 
b quark, and where there may be other jets present (up to two additional jets at NLO). The 
relative importance of Processes 3-6 is significantly increased compared with the exclusive 
cross sections, especially at the LHC, due in particular to Process 6. This is expected 
since Process 6 is more effectively cut in the exclusive cross section where the NLO cross 
section is required to have the same number of jets as the LO cross section. The fixed order 
cross section (Processes 1+2+6) dominates at both the Tevatron and LHC, but corrections 
coming from Processes 3+4+5 are large and of the order of 50% at the LHC. Moreover, 
given the relevance of Process 6, having included part of the NLO corrections to Process 6 
by calculating Process 3 at NLO in QCD increases the stability and therefore the validity of 
the theoretical prediction. The NLO cross sections are also increased by a larger factor than 
in the exclusive cross sections. The fraction of W + Ij + X events in which the jet contains 
a b quark is around 0.9% at the Tevatron and 1.2% at the LHC. 



^As explained in Section [21 Process 6 is the same in both the resummed and fixed order calculations, i.e. 
it does not include any counterterm to subtract the logarithmic terms that have been resummed in Process 
3. Such counterterm has been included with Process 3 since it shares the same final state phase space. 
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Table 2: Exclusive cross sections (pb) for W boson plus one jet, which contains at least one b 
quark, at the Tevatron {^/s = 1.96 TeV pp) and LHC (a/s = 14 TeV pp). No branching ratios 
or tagging efficiencies are included. The labels on the columns have the following meaning: 
Wb = exactly one jet, which contains a b quark; W{bb) = exactly one jet, which contains two 
b quarks; Wj = exactly one jet, including both light quarks (m, d, s, c) and gluons. Both the 
leading-order (in square brackets) and next-to-leading-order cross sections are given. The 
first number given is from Processes 1-2; the second number is from Processes 3-6 (with 
Process 6 given separately in parenthesis); and the third number is their sum. 





Exclusive cross sections (pb) 


Collider 


Wb 


W{bb) 


TeV W+{= W-) 
LHC W+ 
LHC W~ 


[5.28+0.75=6.03] 8.02+0.62(-0.05)=8.64 
[30.2+54.3=84.5] 40.0+48.4(22.6)=88.4 
[21.6+31.4=53.0] 29.8+29.4(12.6)=59.2 


[2.66] 3.73-0.02=3.71 
[17.6] 22.7+11.7=34.4 
[12.9] 17.2+6.5=23.7 




Wj 


TeV W+{= W-) 
LHC W+ 
LHC W- 


[1410] 1790 
[14240] 15810 
[11040] 12040 



We estimate the theoretical uncertainty in the exclusive and inclusive LO and NLO cross 
sections by varying the renormalization (fiR) and factorization (fip) scales independently. 
We vary fi^ in the (0.5 — 2)Mw range while keeping fip fixed at its central value, fip = Mw- 
Similarly we vary fip in the {0.5 — 2)Mw range keeping fiR = Mw- The results are illustrated 
in Figs, m and [5], where the plots on the l.h.s. correspond to W production with a b jet 
(Wb) while the plots on the r.h.s. correspond to W production with a double-6 jet {W{bb)). 
The upper plots refer to W~^b/W~^{bb) = W~b/W~{bb) production at the Tevatron, the 
middle plots to W~^b/W^{bb) production at the LHC, and the lower plots to W~b/W~{bb) 
production at the LHC. The cross sections in Figs. H] and [5] have been normalized to their 
A*R = fJ'F = Myy value. The horizontal axis represents the variation of either fiR or fip, 
depending on the curve (see figure captions), normalized to the central value /xq = M^y. 
In Tables S] and O we quantitatively give the variation with fi^ and /ii? as (asymmetric) 
uncertainties on the central value, corresponding to the choice fiR = fJ^F = Mw used to 
obtain the results of Tables [2] and [31 We have not included in the theoretical uncertainties 
reported in this paper the uncertainty coming from the parton distribution functions. 

From Figs, m (and Eland Tables IH and Ej) we see that the theoretical uncertainty due to the 
dependence on the renormalization scale is larger than the corresponding uncertainty from 
the factorization-scale dependence. The decrease in the factorization-scale dependence in 
going from LO to NLO is mild, while the decrease in the renormalization-scale dependence 
is significant. An exception is the inclusive W{bb) cross-sections at the LHC, where the 
renormalization-scale dependence slightly increases at NLO. Even in the exclusive case, the 
improvement in the renormalization-scale dependence in going from LO to NLO for W{bb) 
is mild at the LHC. These exceptions can be explained by the fact that only Processes 1-2,6 
contribute to W{bb) production at NLO and, among those. Process 6 opens a new initial 
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Table 3: Inclusive cross sections (pb) for W boson plus one jet, which contains at least one 
b quark, at the Tevatron {^/s = 1.96 TeV pp) and LHC {y/s = 14 TeV pp). No branching 
ratios or tagging efficiencies are included. The labels on the columns have the following 
meaning: Wb + X = one or more jets, at least one of which contains a b quark; W{bb) + X = 
one or more jets, one of which contains two b quarks; Wj + X = one or more jets, including 
both light quarks (m, d, s, c) and gluons. Both the leading-order (in square brackets) and 
next-to- leading-order cross sections are given. The ffist number given is from Processes 1-2; 
the second number is from Processes 3-6 red (with Process 6 given separately in parenthesis); 
and the third number is their sum. 





Inclusive cross sections ( 


3b) 


Collider 


Wb + X 


W{bb) + X 


TeV W+{= W-) 
LHC W+ 
LHC W- 


[7.56+1.81=9.37] 11. 77+2.40(0. 77) = 14.17 
[39.3+106.0=145.3] 53.6+136. 1(68. 9) = 189. 7 
[27.9+67.0=94.9] 39.3+88.2(44.6) = 127.5 


[2.66] 4.17+0.39=4.56 
[17.6] 25.1+35.9=61.0 
[12.9] 18.9+23.6=42.5 




Wj+X 


TeV W+{= W~) 
LHC W+ 
LHC W- 


[1410] 2030 
[14240] 20000 
[11040] 15220 



state, namely qg, and is effectively a LO process. The effect is larger at the LHC because, 
due to the large gluon density, the qg channel is more relevant. The effect is also larger for 
inclusive rather than exclusive cross sections because the exclusive final state suppresses the 
contribution of Process 6 (which is a 2 — > 4 process), as evidenced by the numerical results 
in Tables [2] and O 

In Figs. [6H9] we show the differential cross sections with respect to the transverse mo- 
mentum of the b jet and of the W boson, for W~^b inclusive/exclusive and W~^{bb) inclu- 
sive/exclusive production. If there is more than one b jet in the final state, the Pt of the 
highest-pT b jet is chosen. We do not show distributions for W~b and W~{bb) production 
at the LHC, since they resemble the ones for W^b and W~^{bb) production at the LHC il- 
lustrated here. For all final states, the NLO QCD corrections modify the shape of both the 
6-jet and VT-boson transverse momentum distributions. 

4 Conclusions 

In this paper we report on a NLO calculation of the production of a W boson with one b 
jet. We present both inclusive and exclusive cross sections, as well as cross sections where 
the jet contains one or two b quarks. We show that it is essential to keep the fe-quark mass 
finite throughout the calculation, and we are able to overcome this technical hurdle. The 
calculation is performed by combining two previous NLO calculations of Wbb [21] and Wbj 
[22], taking care to treat their overlap consistently. The calculation that we present thus 
represents the state of the art prediction for W plus one b jet production at NLO in QCD. 
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Table 4: Exclusive cross sections (pb) for W boson plus one jet, which contains at least 
one b quark, with the theoretical uncertainty due to renormalization (first uncertainty) and 
factorization (second uncertainty) scale dependence. The uncertainty due to the renormal- 
ization scale (/i/j) dependence is estimated by varying fiji by a factor of two with respect 
to its central value /ir = Mw, while keeping the factorization scale (fip) fixed at its central 
value fip = Mw- The uncertainty due to the factorization scale is estimated analogously, 
i.e. varying fip factor of two about its central value fip = Mw, while keeping fixed 
at fiR = Mw- The central values are extracted from Table O The labeling of columns and 
rows are as described in the caption of Table El Results within brackets are LO, results with 
no brackets are NLO. 



Exclusive cross sections (pb) 



Collider 



Wb 



hU.27+U.U2 

-0.19-0.03 
4-0.27 



)\ 5.b4X(^i„o^;^2-0.003j 



[2.66x(l^-™-) 
[17.6x(ll™;^) 
[12.9x(ll°:?9^1H?) 



W{bb) 



1 0.12+0.01 \ 
-L-O.ll-O.OlJ 
1 +0.23+0.03N 
-■--0.16-0.04 J 
1 +0.21+0.03\ 
-■--0.15-0.04 J 



TeV W+{= W- 
LHC W+ 
LHC W- 



[6.03x(l: 



[84.5x{lt'o-iltl\l)] 88.4x(l+o-i™8 

[53.0X (11:0:^2)] 59^2x(i: 



o.ii-o.io; 

-0.12+0.08 
-0.11-0.10 



3.71x( 
34.4x( 
23.7x( 



These calculations can be compared with the large amount of data on W plus b jets 
already gathered at the Tevatron, and soon to be produced at the LHC. They can also be 
compared with the inclusive samples obtained by merging matrix-element calculations with 
parton showers. 
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Figure 4: Renormalization- and factorization-scale dependence of Wb (l.h.s.) and W{bb) 
(r.h.s.) exclusive production. The upper plots correspond to both W~^b/W^{bb) and 
W~b/W~{bb) production at the Tevatron. The middle plots correspond to W~^b/W^{bb) 
production at the LHC and the lower plots correspond to W~b/W~{bb) production at the 
LHC The black dashed (LO) and red solid (NLO) curves represent the dependence on the 
renormalization scale (fiR) when fiji is varied with respect to its central value fio = Mw, 
while the factorization scale (hf) is fixed aX fip = Mw- In a similar way, the magenta dotted 
(LO) and blue dot-dashed (NLO) curves represent the dependence on np when np is varied 
with respect to its central value /iq = Mjy, while [ipi is fixed at hr = M^. The cross sections 
have been normalized to their fiji = fip = My/ value. 
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Figure 5: Renormalization- and factorization-scale dependence of Wb (l.h.s.) and W{bb) 
(r.h.s.) inclusive production. The upper plots correspond to both W^b/W~^{bb) and 
W~b/W~{bb) production at the Tevatron. The middle plots correspond to W^b/W^{bb) 
production at the LHC and the lower plots correspond to W~b/W~{bb) production at the 
LHC. The black dashed (LO) and red solid (NLO) curves represent the dependence on the 
renormalization scale (fiR) when fiji is varied with respect to its central value fio = Mw, 
while the factorization scale (hf) is fixed aX fip = Mw- In a similar way, the magenta dotted 
(LO) and blue dot-dashed (NLO) curves represent the dependence on np when np is varied 
with respect to its central value /iq = Mjy, while [ipi is fixed at hr = M^. The cross sections 
have been normalized to their fiji = fip = My/ value. 
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Table 5: Inclusive cross sections (pb) for W boson plus one jet, which contains at least 
one b quark, with the theoretical uncertainty due to renormalization (first uncertainty) and 
factorization (second uncertainty) scale dependence. The uncertainty due to the renormal- 
ization scale (/i/j) dependence is estimated by varying fiji by a factor of two with respect 
to its central value /ir = Mw, while keeping the factorization scale (fip) fixed at its central 
value hf = Mw- The uncertainty due to the factorization scale is estimated analogously, 
i.e. varying fip factor of two about its central value fip = Mw, while keeping fixed 
at fiR = Mw- The central values are extracted from Table [3l The labeling of columns and 
rows are as described in the caption of Table [3l Results within brackets are LO, results with 
no brackets are NLO. 





Inclusive cross sections (pb) 


Collider 


Wb 


W{bb) 


TeV W+{= W-) 
LHC W+ 
LHC W- 


[9.37x(ll™:[{^)] 14.17x(lt[;;K^r) 
[U5.3xilt'oiltrd)] 189.7x(nin?o^) 
[94.9x(lt0:?^t°;}^)] 127.5x(ll°:in?^) 


[2.66x(ll-™)] 4.56x(l^-™) 
[17.6x(lt™?^)] 6L0x(lt°ft°:0^) 
[12.9x(lt™??)] 42.5x(lt°ft°:°i) 
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Figure 6: prib) (upper plots) and priW) (lower plots) distributions for W~^b exclusive pro- 
duction, at the Tevatron (l.h.s. plots) and at the LHC (r.h.s. plots). 
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Figure 7: prib) (upper plots) and priW) (lower plots) distributions for W~^{bb) exclusive 
production, at the Tevatron (l.h.s. plots) and at the LHC (r.h.s. plots). 
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Figure 8: prib) (upper plots) and Pt{W) (lower plots) distributions for W~^b inclusive pro- 
duction, at the Tevatron (l.h.s. plots) and at the LHC (r.h.s. plots). 
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Figure 9: prib) (upper plots) and priW) (lower plots) distributions for W^{bb) inclusive 
production, at the Tevatron (l.h.s. plots) and at the LHC (r.h.s. plots). 
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